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ABSTRACT 

We explore how well crowded field point-source photometry can be accomplished with Sloan Digital 
Sky Survey (SPSS) d ata. For this purpose, we present a photometric pipeline based on DoPHOT 
(jSchechter et al.lll993l ). and tuned for analyzing crowded-field images from the SDSS. Using Monte 
Carlo simulations we show that the completeness of source extraction is above 80% to an i band AB 
magnitude of < 21 and a stellar surface density of ~ 200 arcmin -2 . Hence, a specialized data pipeline 
can be efficiently used for fairly crowded fields, such as nearby resolved galaxies in SDSS images, 
where the standard SDSS photometric package Photo, when applied in normal survey mode, gives 
poor results. 

We apply our pipeline to an area of ~ 3.55D° around the dwarf spheroidal galaxy (dSph) Leo I. 
Using the resulting multi-band (g,°,i) photometry we construct a high signal-to-noise star-count map 
of Leo I, utilizing an optimized filter in color-magnitude space. This filter reduces the foreground 
contamination by ~ 80% and enhances the central stellar surface density contrast of the dwarf by a 
factor of > 4, making this study the deepest wide-field study of the Leo I dSph based on accurate 
CCD photometry. We find that the projected spatial structure of Leo I is ellipsoidal. The best fitting 
empirical King model to the stellar-surface density profile yields core and tidal radii of (6.21 ± 0.95)' 
and (11.70 ± 0.87)', respectively. This corresponds to (460 ± 75) pc and (860 ± 86) pc assuming a 
distance to Leo I of 254j^g kpc. The radial surface-density profile deviates from the King profile 
towards outer radii, yet we find no evidence for 'S' shaped or irregular tidal debris out to a stellar 
surface-density of 4x 10~ 3 of the central value. From the luminosity function of all possible Leo I stars, 
which we carefully extrapolated to faintest magnitudes, we determine the total I c -band luminosity of 
Leo I to be (3.0 ± 0.3)xl0 6 L/ C! ©. We model the mass of the dSph using the spherical and isotropic 
Jeans equation and infer a central mass density of 0.07 M Q pc~ 3 leading to a central mass-to-light ratio 
of ~ 3 in I c band solar units. Assuming that the mass in the system follows the distribution of the 
visible component, we constrain a lower limit on the total mass of the dSph to be (1.7±0.2)xl0 7 M©. 
On the other hand, if the mass in Leo I is dominated by a dark-matter (DM) halo with constant 
density, then the mass within the central 12' yields (2 ± 0.6) xlO 8 M©. Combining the inferred mass 
estimates with the total luminosity leads to a mass-to-light ratio of 3> 6 in I c band solar units, and 
possibly > 75 if the DM halo dominates the mass and extends further out than 12'. In summary, our 
results show that Leo I is a symmetric, relaxed and bound system; this supports the idea that Leo I 
is a dark-matter dominated system. 

Subject headings: surveys; galaxies: Local Group, dwarf, halos; cosmology: dark matter; methods: 
data analysis; techniques: image processing 



1. INTRODUCTION 

1.1. The Local Group as seen by the Sloan Digital Sky 

Survey 



Studies of the Local Group have been recently in- 
vigorated by sensitive multi-wavelength large-area sur- 
veys, in particular the Sloan Digital Sky Survey (SPSS; 
York e t al. 200 0UStoughton et a l. 2002: Abaz aiian et all 
20031 [20041 . [20051 : lAdelman-McCarthv et al.ll2006h . To 
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date SDSS has publicly released high-quality near UV to 
near IR five-ban d photometry and accurate astrometry 
(|Pier et al.! 2003) for ~215 million objects sel ected over 
8000 □ °(DR5: lAdelman-McCarthv et al.ll2007h . The sur- 
vey utilizes highly-automated reduction packages. In 
particular, the pipeline which extracts the multi-band 
photometry from SDSS images is call ed Photo (a detailed - 
descript ion of the pipeline is g iven in iLupton et al]|2002t 
see also iStoughton et"aLll2002f ). 

The SDSS has had a significant impact on Local Group 
studies. It h as led to the discovery of the faintest 
known dSphs ( WillmaneLaLl 120051 : IZucker et aD 120061 : 
iBelokurov et al.1 12006al fbl|d1) as well as to new insights 
into already known dSp h galaxies (e.g. Draco dSph; 
lOdenkirchen et alj 1200 ll) . Analyses of the SDSS pho- 
tometric dataset have revealed tidal streams surround- 
ing t h e Milky Way (for example. lOdenkirchen et al.l 
20011: lYannv et al.l 120031: iGrillmair fc Dionatosl 1200a 
Belokurov et al.1 l2006aU blldl). There have b een detailed 
struct u ral studies of the Galacti c halo (|Helmi et al.l 
2003al IXu et all 120061) and disk (|Helmi et alj l2003bt 
Juritsch et al. based on SDSS. However, there is a 
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serious problem associated with the SDSS standard pho- 
tometric pipeline: Photo does not extract photometry in 
crowded fields such as globular clusters or nearby galax- 
ies (see ij2] for details). For example, the photometry is 
incomplete for the SDSS fields 1 at the centre of the Leo I 
dSph galaxy (see Fig. [T]). 

The aim of this paper is twofold. First, we present a 
photometric pipeline targeting SDSS images of crowded 
stellar fields. This pipeline i s based on DoPHOT, a soft- 
ware package developed by ISchechter et al"] fl993), and 
is designed to be highly automated. Secondly, we apply 
this pipeline to a 3.55D° area centered on the Leo I dSph, 
and study the properties of the dwarf. 

1.2. Leo I dwarf spheroidal galaxy 

The dwarf spheroidal (dSph) ga laxy Leo I was discov- 
ered in the Palomar Sky Survey bv lHarrington fc Wilson! 
(1950). Subsequent studies using photographic plates 
were substantially hampered by the nearby (~20' south 
of the dSph's center) first magnitude foreground star 
Regulus (a Leonis). As a result, the first Leo I 
color-m agnitude diagram (CM P) was published only re- 
cently (|Fox fc Pritchetl 119871) . Nevertheless, Hodge & 
Wright (1978) had already reported the presence of an 
intermediate-age stellar population in the galaxy in- 
dicated by an unusually large number of anomalous 
Cepheids. Later studies based on CCD photometry in- 
dicated that the intermediate age sta rs are the domi- 
nant population dReid fc Mould! Il99lt iLee et all fl99l 
iDemers et~aT1 1 1 9941 ). Base d on Hubble Space Tele- 
scope (HST) observ ations, ICaputo et al.1 (|1999j ) and 
iGallart et al.l (1999a, b) showed that Leo I has an ex- 
tended star formation history. The oldest stars in Leo I 
formed approximately 9-13 Gyr ago, while the youngest 
stars are less than 1 Gyr old. It is this latter trait which 
makes Leo I unusual compared to most other dSphs. 
An old stellar population with an age > 10 Gyr has 
been found just recent ly in the outer regions of Leo I 
dHeld et al.ll2000l l200lh . 

An SDSS field is defined as an imaged area in the sky consisting 
of five (u,g,",i,z) frames. 
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Fig. 1. — Stellar density plots of Leo I shown in Great circle 
coordinates. The grey-scale is logarithmic, and the same in both 
panels; darker colors display higher densities. The top panel shows 
the SDSS (Photo) photometry of Leo I. The bottom panel shows 
DoPHOT's photometry for the same area. Both panels show only 
objects classified as stars by each individual package. Note the 
missing stars in and around the center of the galaxy in the top 
(SDSS - Photo) panel. 

Leo I is thought to be one of the most distant satel- 
lites of the Milky Way i n the Local Group. T he most 
recent distance estimate (jBellazzini et al.ll2004| ) puts the 
galaxy at 254^ig kpc away from our Galaxy [they infer 
a distance modulus of (m — M)o = 22.02 ± 0.13]. 

Generally, the high velocity dispersions for Local 
Group dSphs, combined with dynamical mass estimates, 
indicate that these systems are dominated by dark mat- 
ter (DM; for a review see iMateol Il998al ). The stan- 
dard method for estimating mass-to-light (M/L) ratios 
of pressure-supported systems is to use one-component 
isotropic King models (King 1962, 1966; Richstone & 
Tremaine 1986) with simplifying assumptions that the 
stellar velocity dispersion is isotropic and, critically, that 
mass follows light. Recent kinematic studies have weak- 
ened this latter assumption; velocity dispersion profiles 
do not demonstrate the characteristic decrease with ra- 
dius predicted by isotropic King models. Hence it is 
likely that 'traditional' M/L estimates for dSphs are in 
fact lower than the true value. Specifically, iMateo et al.l 
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(|1998bf) estimate the Fband M/L ratio of Leo I to lie 
in the range of 3.5 - 5.6 (in solar units; see also iMateol 
Il998al) . This result is based on the central velocity dis- 
persion measured for 33 Leo I stars using the structural 
parameters (such as the core and tidal radii) given in 
Irwin & Hatzidimitriou (1995; hereafter IH95). 

IH95 determined the morphology of eight Local Group 
dSphs from star counts using digitized photographic 
plates. One major advantage of this study (compared 
to earlier studies which were largely based on eyeball 
counts) was that an objective star count method was ap- 
plied to all dwarfs. However, the usage of photographic 
plates is for obvious reasons much more restricted then 
the usage of CCD photometry, which we utilize here. 
One general problem, which was not taken into account 
by IH95, is the effect of star blending, which may be 
significant in regions of high number densities such as 
the cores of compact stellar systems. Thus, blending 
may present a serious problem in the determination of 
structural properties of dSphs as the number densities 
are not correctly estimated. Here we carefully approach 
this problem in order to obtain the correct morphological 
parameters of the Leo I dSph. In summary, better con- 
straints of the structural properties of dSph galaxies are 
essential for more robust derivations of masses and M/L 
ratios. In particular, given the large dista nce from the 
Milky Way and the great systemic velocity (jMateo et al.1 
1998b), Leo I plays a cruc ial role in the estimates of the 
mass of our Galaxy (e.g. IZaritskv et al.1 1989: Z aritskvl 
119991 ) and the whole Local Group (|Lvnden-Bel]|ll999h . 

With the advent of multi-wavelength large-area sky 
surveys, such as the SDSS, the properties of nearby 
galaxies can be comp rehensively investigated (e.g. 
lOdenkirchen et ail 2001). In this paper we focus on de- 
riving structural properties of Leo I utilizing SDSS im- 
ages. We use these to constrain more robustly the mass 
and total luminosity of Leo I and finally model the mass 
of the dSph. We also address the total mass estimates 
for Leo I with and without the mass-follows-light assump- 
tion. Finally, we derive a range of possible M/L ratios of 
the galaxy utilizing the improved structural parameters. 

Our photometric pipeline is described in Sj2] The exten- 
sive tests performed on the pipeline in order to infer the 
photometric accuracy and completeness are presented in 
The results on Leo I are given in §^3HSJ SI describes 
the color-magnitude selection of Leo I candidates, i.e. the 
construction of a high-contrast map of Leo I, and briefly 
the color-magnitude diagram (CMD) of the dSph. We 
constrain the size and structure of Leo I in 351 an( l de- 
rive the total luminosity, mass and M/L ratio in $6l We 
discuss our results in [J7] and summarize them in EjSj 

2. THE PHOTOMETRIC PIPELINE 

The standard SDSS photometric pipeline, Photo, 
has two main restrictions: (i) it is time limited (to 
1 ms/object), and (ii) the number of extracted objects 
per image cannot exceed a given number. The photome- 
try package therefore does not analyze centers (and some- 
times whole frames) of crowded regions (for example, 
nearby resolved galaxies). This is illustrated in Fig. [TJ 
Here we present a pipeline developed for obtaining high 
quality photometry from SDSS images of crowded fields. 
The pipe line is based on the DoPHOT source extraction 
package (jSchechter et al.lli"993l ) designed to search for ob- 



jects on a digital image of the sky and to extract posi- 
tions, magnitudes and classifications for those objects. 
The package is widely used in the astronomical commu- 
nity and has a pr o ven accuracy (e.g.lReid fc Mouldlll991l: 
IVoet et all 119951 : iGallart et all Il999bl : iBellazzini et all 
12004 ). We use a version of DoPHOT 2.0 which was 
slightly modified by Eugene Magnier. 2 

In the following work we focus on g, ° and i bands since 
the u and z bands are less sensitive. Our pipeline can 
be divided into four parts: 1) aligning the SDSS frames 
and extracting coefficients (from the SDSS tsField 3 file) 
needed for further calculations, 2) adjusting DoPHOT's 
input parameters, 3) running DoPHOT on the individual 
frames and 4) converting the output photometry to the 
AB magnitude scale and computing several other quan- 
tities, as well as extracting the proper astrometry for a 
given field (see £12.21 for details). We briefly describe the 
pipeline below. 



2.1. DoPHOT's input parameters 

DoPHOT requires an initial list of parameters (fine- 
tuned for the given image) to run. In order to make 
the object-extraction from SDSS images fast and highly 
automated we developed template DoPHOT parameter 
files for each photometric band. We fine-tuned these pa- 
rameters and finally adopted the ones giving the best 
detection rates. 

Two sets of parameters required close attention: the 
'aperture box' values, and those defining the background 
sky model. Since DoPHOT fits a model point spread 
function (PSF) to the data (as opposed to a numerically 
given empirical PSF), it also computes aperture fluxes to 
correct for the systematic errors introduced by using a 
model PSF. These 'aperture box' parameters were fine- 
tuned to reproduce the existing SDSS photometry in un- 
crowded fields, where Photo performs well 4 . The second 
important parameter set was those defining the back- 
ground sky model. Crowded fields result in a background 
sky level which varies significantly across the field, and 
DoPHOT can model this variation using either a uniform 
gradient model or a modified Hubble profile (for the lat- 
ter see e.g. Binney & Tremaine 1987). We adopted the 
Hubble model since it increased the detection rate by a 
factor of ~3 due to the more realistic background sky 
description. 

These two sets of parameters were then hard- wired into 
the fine-tuned template input list. Yet, there are two val- 
ues that need to be specified separately for each frame, 
namely estimates of seeing and background sky. We cre- 
ated the pipeline in such a way that it measures these 
two values directly from the image and adds them to the 
parameter template to create the final input files. 

2 The entire program is wrapped inside a C program which 
implements dynamic memory allocation, furthermore the C code 
interprets command-line arguments and allows compilation under 
f2c. 

3 An SDSS tsField file is a binary FITS table which contains pa- 
rameters relevant for the entire field (for details see www.sdss.org). 

4 The parameters were tested on fields surrounding the Leo I 
dSph. Nonetheless, if the need for different aperture corrections 
(or any of the fixed parameters) arises, the parameters can easily 
be changed and added to the pipeline without adversely affecting 
its performance. 
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2.2. Calibrating DoPHOT's output to the AB and 
J2000 systems 

During its photometry routine, DoPHOT outputs a list 
of detected objects with several quantities, including the 
object positions and classification, total and 'aperture' 
magnitudes (with corresponding errors). DoPHOT's 
magnitudes are given in the form of —2.5 log (DN), where 
DN are counts in units of digital numbers. To convert 
them to the AB 5 magnitud e scale (for details ab out the 
AB photometric system see lFukugita et al1l!996l and ref- 
erences therein) we used the parameters from the SDSS 
tsField files to calibrate DoPHOT counts: 



-2.5 log 



\t c 



CO - C oxt • TOair 



(1) 



where i eX p is the exposure time (t cxp = 53.907456 s), |co| 
is the magnitude that yields DN/i exp of 1 at zero airmass, 
Cext the extinction coefficient and m a i r the airmass. The 
values co, c cx t, rn^ with their uncertainties are extracted 
from the SDSS tsField file 6 . 

A final magnitude correction is calculated via a set 
of 'perfect' stars, those objects which best represent 
the PSF shape. [To find these stars DoPHOT iterates 
through the given image starting at the highest level. 
At a given level, it creates a typical stellar shape using 
all detected stars within this iteration. Then it fits this 
shape to a given object and determines the magnitude. 
This procedure is repeated at each level. If the shape and 
magnitude fits converge for a given object, it is classified 
as a 'perfect' star.] The final correction is then calculated 
as the median offset between the total and aperture mag- 
nitude systems for these perfect stars. The correction is 
applied to the complete object list, providing a final set 
of total AB magnitudes (calculated from the total fluxes 
derived from fitting a model PSF). 

The AB magnitude errors are calculated taking into 
account the errors of the variables used in eq. [T] . Those 
are then added to the median uncertainty in quadrature 
to calculate the final uncertainties. Thus, for any given 
object reported by DoPHOT, the magnitudes are scaled 
to the AB system, aperture corrected and the uncertain- 
ties are reported. This is completed separately for each 
frame. From the positions in pixel coordinates as given 
by DoPHOT the J2000 right ascension and declination 
for each object are recovered. The objects extracted in 
each frame are then matched between all bands allowing 
for a maximum separation of ~ 2.8". Thus, a final stellar 
list with photometric magnitudes and uncertainties in g, 

and i is obtained. In order to compare this photome- 
try with the existing SDSS data, the match algorithm is 
used to find common stars between the two datasets. 

3. TESTING DoPHOT 

In this section we describe tests performed on SDSS im- 
ages using our automatic package (see previous section). 
Our aim was to constrain the photometric accuracy and 

SDSS photometry is intended to be on the AB system 
(Okc & Gunn 1983; Fukugita ct al. 1996), by which a magnitude 
object should have the same counts as a source of F v = 3631 Jy. 
However, the photometric zeropoints are slightly off the AB stan- 
dard. Nonetheless, according to the present estimate the (g,°,i) 
band zeropoints are close to the AB system (~ 0.01 mag). For 
details see www.sdss.org. 

6 The SDSS tsField files report c < 0. 



the detection efficiency of the pipeline, in particular for 
crowded fields where the SDSS photometric software, 
Photo, gives poor results (see Fig. [TJ. The results below 
demonstrate that our DoPHOT-based pipeline can ex- 
tract high quality photometry from SDSS crowded-field 
images. 

3.1. Photometry 

The accuracy of DoPHOT's (g, 5 i) magnitudes rela- 
tive to SDSS imaging magnitudes is summarized in Fig. [2] 
for stars in an un-crowded field. The left column shows 
the magnitude difference as a function of SDSS magni- 
tude. The middle column shows the histogram of the 
magnitude difference at the bright end (mag < 21). The 
distribution of the magnitude difference is symmetric in 
the ° band, while the g and i bands show slight sys- 
tematic trends caused by the differences between the 
two pipelines, Photo and DoPHOT. One possible source 
for this skewed behavior may be the aperture param- 
eters which were used as an input for DoPHOT, and 
were optimized in the r band. It is also possible that 
the different sky background estimations used in the two 
pipelines contribute to these effects. The right column 
in Fig. [2] displays the magnitude differences normalized 
by the expected errors determined by adding the SDSS 
and DoPHOT errors in quadrature. Note, that the er- 
rors include both photon noise errors and other system- 
atic uncertainties. The mean equivalent Gaussian width, 
<7, indicated in each panel, is close to unity and demon- 
strates that various possible errors in extracting counts 
do not exceed the noise expected from photon statistics. 

3.2. Detection efficiency 

In order to estimate the accuracy of DoPHOT's pho- 
tometry and the efficiency of recovering stars as a func- 
tion of both magnitude and stellar density, we performed 
artificial star tests in crowded and un-crowded fields. 
The tests were performed separately for (g, ; i) bands. 
Synthetic stars were generated using the DoPHOT PSF 
model, which consists of similar ellipses of the form: 



I(z)=h 
1 



3 1 



k\ 



.k=0 



(2) 
(3) 



To avoid over-crowding (i.e. blending), the artificial 
star tests were done via Monte Carlo simulations by 
adding 100 randomly distributed synthetic stars with an 
appropriate FWHM into each frame. To assess the pho- 
tometric accuracy of the recovered stars, the procedure 
was repeated 100 times. The photometry was then ex- 
tracted utilizing our pipeline as described in Sec. [2j In 
Fig. [3] we summarize the photometric accuracy for the 
same un-crowded field shown in Fig. O The results for 
a crowded field are shown in Fig. |4j The magnitudes 
are reproduced with the expected accuracy (a w 0.014), 
with no systematic errors in the overall magnitude scale 
(the width of the distribution in the last column is very 
close to unity). The slight systematic offset in the mag- 
nitude differences of « 5xl0 -4 and w 5xl0 -3 for the 
un-crowded and crowded field, respectively, is due to the 
aperture correction (the 'aperture box' parameters were 
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Fig. 2. — The comparison of DoPHOT and SDSS magnitudes 
in an un-crowded field (run 3631, camcol 3, rerun 40, field 236) in 
the surrounding of Leo I. The first column displays the magnitude 
difference as a function of magnitude; the large symbols connected 
by lines show the ±3<r envelope. The distribution of magnitude dif- 
ferences at the bright end is displayed in the middle column. The 
last column shows the distribution of the magnitude differences 
normalized by the expected errors of DoPHOT's and Photo's pho- 
tometry added in quadrature at the bright end (solid line) and for 
the full sample (dashed line). 

set to fit SDSS magnitudes in un-crowded regions). Since 
the offset is of order of 10~ 3 mag, while the magnitudes 
are accurate to 10~ 2 mag, the discrepancy is negligible. 
We conclude that our photometry is sufficiently accu- 
rate for robust studies of, for example, nearby galaxies 
resolved by SDSS. 

For a thorough analysis of star counts, luminosity func- 
tions, surface density profiles etc., one needs to under- 
stand the limits of the source extraction software. It 
is already well known that detection efficiency declines 
with stellar brightness. However, since this work is fo- 
cused on crowded fields, it is essential to understand how 
the number of retrieved stars depends on the local stel- 
lar density. To quantify these dependences we define the 
completeness as the ratio of the number of artificial stars 
extracted by DoPHOT to the number of artificial stars 
added to the frame, n ou t pu t /"-input ■ Fig.[5]shows the com- 
pleteness as a function of magnitude for different SDSS 
fields. We consider a star as recovered if i) the position 
as extracted by DoPHOT matches the position of the 
inserted artifical star (within a box of 2" on a side cen- 
tered on the star) and ii) the classification for that object 
given by DoPHOT is a 'perfect' star. The uncertainties 
in Fig. [5] were derived from Poisson errors for the input 
and output star counts in 1 mag wide bins. Therefore, 
they reflect the uncertainties due to the low number of 
artificial stars injected rather than the intrinsic uncer- 
tainty of the completeness as a function of magnitude. 
The fraction of recovered stars falls below 90% at mag- 
nitudes fainter than 20 — 21 (depending on the filter). 
For reference, the quoted SDSS 95% completeness is in 
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Fig. 3. — Artificial star test results for an un-crowded field 
(same field as in Fig. [2j. The first column displays the magnitude 
difference of injected and extracted magnitudes as a function of 
magnitude; the line connected with symbols represents the ±3cr 
envelope. The middle column shows the distribution of magnitude 
differences for mag< 21. The last column displays the distribution 
of magnitude differences normalized by the expected errors at the 
bright end (solid line) and for the full sample (dashed line). 
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Fig. 4. — Same as Fig. [4] but for a crowded field encompassing 
a part of Leo 1 (run 3631, camcol 3, rerun 40, field 238). 

the range of 21.3 — 22.2 mag for the (g,r,i) bands (using 
Photo); note, however, that this value is only represen- 
tative for un-crowded regions. 

Our next aim was to measure the completeness as a 
function of stellar density. The following tests were per- 
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Fig. 5. — The fraction of recovered artificial stars (completeness) 
as a function of magnitude for an un-crowded field surrounding 
Leo I (see text for details). Note that the completeness starts 
declining rapidly for mag > 21. 

formed on two crowded fields at the center of the Leo I 
dSph (run 3631, camcol 3, rerun 40, fields 237 and 238). 
Each field was divided into cells of 50 x 50 square pix- 
els in area, and the stellar density in each cell was then 
measured prior to artificial star injection. The artificial 
stars were then placed in the field and recovered using the 
pipeline, and the completeness fraction was measured for 
each cell. To ensure we obtained accurate completeness 
measurements, the Monte Carlo simulations described 
above were executed for 1000 iterations. Thus, 10 5 stars 
were randomly placed in each individual (g, ; i) frame. 
The combined results are displayed in Fig. [6] Note that 
we show the completeness for four magnitude bins. The 
crowding at the center of Leo I has no significant effect 
on stellar recoverability for stars brighter then ^20th 
magnitude. Stellar crowding reduces the completeness 
at fainter magnitude bins, inflicting a 10 — 30% loss for 
stellar densities up to ^200 stars/arcmin 2 . 
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Fig. 6. — The fraction of recovered artificial stars (completeness) 
as a function of stellar density for different magnitude bins. The 
panels present the combined artificial star tests ran on fields 237 
and 238 (run 3631, camcol 3, rerun 40). The band for which we 
show the derived completeness is indicated in each panel, and each 
curve represents the completeness for one magnitude bin, which 
is indicated in the legend on the righ-hand side. Note that the 
completeness is essentially not a function of stellar density till 
mag ~ 20. 

In Table [1] we summarize the completeness as a func- 
tion of both magnitude and stellar density for the i band. 
In latter sections we will use these results to correct star- 
counts and flux integrals for incompleteness introduced 
by DoPHOT. 

4. PHOTOMETRIC PROPERTIES OF Leo I 

In this section we present the construction of a high- 
contrast map of Leo I utilizing the color-magnitude in- 
formation, briefly describe the main features of the color- 
magnitude diagram (CMD) of Leo I, and compare them 
to previous results. 

4.1. Photometric selection of Leo I candidates 



Table [2] summarizes the SDSS great circle scans which 
encompass ~ 3.55D around the Leo I dSph galaxy (see 
Fig. [?])■ The photometry extraction was performed as 
descri bed in [J2] An examination of the Schlegel ct al. 
(1998) dust map shows that extinction is not a signifi- 
cant problem for this object; reddening varies by approx- 
imately 0.01 mag across the Leo I field. The mean value 
at the centre of the dSph is E{B — V) = 0.036, and thus 
the entire photometric dataset was corrected using this 
value. The empty area south of Leo I in Fig. [7] is a re- 
sult of contamination by the 1st magnitude star Regulus. 
This and other regions around bright stars are excluded 
in the further analysis. 
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Fig. 7. — Spatial distribution of objects classified as 'perfect' 
stars by DoPHOT. The shown area covers ~ 3.55d° around Leo I. 
The dSph is visible in the central region of the panel. The top panel 
shows the distribution in right ascension and declination (J2000), 
the bottom in great circle coordinates (used for scanning the sky 
by the SDSS). The empty area in both panels, south of Leo I, is due 
to contamination by the first magnitude foreground star Regulus 
(a Leonis). 

In order to determine the photometric characteristics 
of the Leo I population, we selected stars within an ellipse 
of 6' semi-major axis (approximately half the limiting ra- 
dius) in the central region of Leo I. The distribution of 
field stars was obtained from stars beyond an ellipse with 
a semi-major axis of 30'. The color-magnitude diagram 
(CMD) for stars in the central region of Leo I is shown in 
the top left panel in Fig. [5[ where the field star distribu- 
tion is overlaid in contours. As anticipated, the distribu- 
tion of field stars shows two main populations: the bluer 
halo stars (g — r f» 0.4) and the redder (g — r « 1.35) 
thick disk stars. 
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TABLE 1 

i BAND COMPLETENESS AS A FUNCTION OF MAGNITUDE AND STELLAR DENSITY 



density I fraction [%] upper lower fraction [%] upper lower fraction [%] upper lower fraction [%] upper lower 
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TABLE 2 
Specifications for SDSS images 

UTILIZED IN THIS PAPER 



Run 


Rerun 


Camcol 
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To robustly study the dSph galaxy, we needed to con- 
struct a map which maximized the number of Leo I 
stars and minimized the level of field star contamina- 
tion. In order to do this, we followed the method of 
iGrillmair et alJ (119951). which was re fined for SDSS pho- 
tometry bv lOdenkirchen et al.l ((2001) in their analysis of 
the Draco dSph. Stars within the inner region of Leo I 
form a tight correlation in the g — r vs. r — i plane. This 
is illustrated in the top right panel in Fig. [7] To quantify 
the width of this distribution we defined a set of prin- 
cipal axes (ci, C2) where c\ measures the position along 
the g — r vs. r — i locus, and C2 the position perpendicular 
to it: 

ci = 0.920 • (g - r) + 0.391 ■ (r - i) (4) 
c 2 = -0.391 • (g - r) + 0.920 • (r - i) (5) 

The bottom panels in Fig. [8] show the CMDs for stars 
within the inner Leo I region using the principal colors. 
Since c% is essentially a measure of the photometric dis- 
persion, we imposed a cut of C2 ^ 2er C2 to Leo I stars. 
Here a C2 designates the standard deviation of the C2 dis- 
tribution measured as a function of i magnitude (see bot- 
tom right panel in Fig. [5]). We imposed an additional 
magnitude cut of i < 22 to the stars in the analyzed area 
since the completeness decreases significantly at fainter 
magnitudes and the photometric uncertainties are much 
higher. 

Following the method of lOdenkirchen et alJ (|2001l ). we 
constructed CMD functions for the field and Leo I pop- 
ulations. The (ci,i) CMD was divided into a series of 
cells, where the cell dimensions were 0.09 mag in color 
and 0.35 mag in z-magnitude. Additionally, each cell was 
separated from its 'neighbors' by 0.015 mag in color and 
0.05 mag in i. This level of overlap between neighboring 
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Fig. 8. — Color-magnitude and color-color diagrams of stars in 
the region around Leo I. In the top left panel dots represent the 
central region of Leo I, i.e. stars within an ellipse with a semi- 
major axis of 6' in the central region of Leo I. Overlaid contours 
show the distribution of field stars obtained from stars beyond an 
ellipse with a semi-major axis of 30' centered at Leo I (see text 
for details). The magnitudes are corrected for reddening. The 
top right panel shows the r — i vs. g — r diagram for stars with 
i > 21, which were used to obtain the principal component axes 
(ci, C2), which are indicated in the panel (see eqs. 0] & [5] and 
text for details). Note that the position of the crossing point of the 
principal axes was shifted for clarity. The bottom panels show the 
color-magnitude diagrams for the stars in the inner Leo I region 
(dots) utilizing the principal component colors. The contour in the 
bottom left panel indicates the optimized photometric filter for the 
selection of potential Leo I members. The connected large points 
in the bottom right panel indicate the ±2<r envelope of the C2 vs. 
i distribution. 

cells ensured the CMD function was smooth and contin- 
uous. As described above, the Leo I stellar population 
was selected from the inner 6' of the dwarf, while the field 
population was taken from the area beyond a radius of 
30'. This second value is ~2.5 times the Leo I tidal ra- 
dius, and hence the field area is not expected to contain 
any of the dSph's stars. The CMD function was then 
constructed for both the field and Leo I po pulati ons. 
As described by IGrillmair et al.l (|1995f) and 
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lOdenkirchen et al] (|200lD , we then derived a 'sig- 
nal' for each CMD cell by comparing the Leo I and field 
populations. For a given cell at position (i,j), the signal 
was calculated as: 



n c (i,j) -gn t (i,j) 



(6) 



where n c (i,j) and nt{i,j) describe the Leo I (core) and 
field populations respectively. The factor g is a scaling 
factor defined as the ratio of the core to field areas. We 
optimized the population contrast s using a threshold 
value So such t hat s ^ Sp. The va l ue of So was derived 
using eq. [2] in lOdenkirchen et all (|2001h in such a way 
that the contrast was optimized for an annulus between 
30' and 60' from the center of Leo I. The optimal color- 
magnitude filter mask defined by so is outlined in the 
bottom left panel in Fig. [5J This filter removed ~ 80% 
of the field stars and enhanced the central stellar density 
contrast by a factor of > 4. 

4.2. Color-magnitude diagram 

The morphology of the CMD presented here (see left 
panels in Fig. [8j agre es very well with those obtained by 
previous studies (e.g. [Reid fe Mouldll99lllDemers et alJ 
fl99ilGallart et al.lfl999al lbl: lHeld et al.ll2000ll200lD . The 
most prominent feature is a well defined red giant branch. 
The red clump is visible in the range 0.1 < g — r < 1 and 
21 < i < 22.5. The stars in the horizontally extended 
region at i « 20.5 and —0.3 < g — r < 0.3 correspond to 
anoma lous Cepheids reported first by Wallerstein fe Coxl 
and explored further by e.g. iLee et al.l (|1993l ). 



5. SIZE AND STRUCTURE OF Leo I 
5.1. The center 

In this section our aim is to quantify the size and struc- 
ture of Leo I. Hence it is important to constrain the 
center of the stellar surface-density distribution of the 
dSph. In order to do so we spatially binned the color- 
magnitude filtered Leo I candidate stars and fitted a 2- 
dimensional Gaussian to this distribution. To account for 
any artifacts due to our procedure and to robustly esti- 
mate the uncertainties, we derived the center, position 
angle (PA) and ellipticity 7 of the dwarf using bootstrap- 
resampling. The mean center in great circle coordinates 
and the corresponding errors are /i = (152.970 ± 0.003)° 
and v = (-0.365 ± 0.002)°. The best fit position angle 
for great circle coordinates is PA = (— 1.7 ± 5.7)° and 
the ellipticity is 0.3 ± 0.1. In the remainder of the paper 
we will take these coordinates as the center of Leo I and 
use the above derived PA and ellipticity values to deduce 
the surface-density profile of the dSph. 

The right ascension and inclination of the ascending 
node of the great circle scans of the area analyzed here 
are 95° and 15°, respectively. Hence, the center in 
J2000 coordinates is at a = 152.122° ± 0.003° and 5 = 
12.313°±0.003°, i.e. 10 ,t 08 m 29.4 s ±0.8 s + 12 d 18'48"±9", 
the PA is (-9.2 ± 5.7)° and the ellipticity stays the 
same. Our PA and ellipticity agree within the uncertain- 
ties with the values given in e.g. IH95. The center de- 
rived here is roughly consistent with the values from past 

7 Ellipticity is computed as 1 — b/a, where a and b are the major 
and minor semi-axes of the ellipse. 



studies (IZwickv et al.lll961t iGallouet fe Heidmanr] |1971l: 
iDressel fc Condon! 119761 : iFalco et all 119991 : see NED for 
a summary). It is worth noting, however, that the posi- 
tions of the center of Leo I vary in the literature by ~ 20' 
and ~ 30' in right ascension and declination, respectively. 

5.2. The size 

In Fig. [9] we show the projected density profile of the 
CMD-selected stars. The radial stellar surface-density 
was obtained by computing stellar densities within ellip- 
tical annuli of 1.2' width starting at the center of Leo I 
(the center, ellipticity and position angle were chosen 
as described above). Within each annulus we corrected 
for incompleteness (see Tab. Q}, thereby accounting for 
the brighter photometric limit of the crowded central re- 
gion (see also Fig. [5]). The datapoints in Fig. [HI display 
the completeness-corrected density (E) above the back- 
ground (also CMD selected and corrected for incomplete- 
ness; Ebkg = (0.340 ± 0.005) arcmin -2 ) as a function of 
the radius (i.e. the rms of the outer and inner semi-major 
axis for a given annulus). The projected background den- 
sity was estimated using the distribution of field stars 
beyond an ellipse of semi-major axis of 30' and the same 
ellipticity as above centered at Leo I. The error bars in- 
dicate quadratically combined S and Ebkg errors, where 
both take into account completeness and Poisson uncer- 
tainties. 
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Fig. 9. — Radial profile of the completeness-corrected and 
background-subtracted stellar density for potential Leo I members 
(dots) . The error bars indicate quadratically combined errors of the 
background and stellar surface-densities, where both are taken to 
be dominated by Poisson and completeness uncertainties. Shown 
are the best-fit King (King 1962; dashed line), Plummer (Plummer 
1911; dash-dotted line), and Gaussian (solid line) profiles. We use 
the latter to estimate the mass of Leo I (see text for details). The 
best fit parameters for all fitted models are indicated in the panel. 
Note that the Plummer profile gives the poorest fit to the data. 



To quantify the size of Leo I we fit an empirical King 
profile (King 1962; hereafter K62) to the stellar surface- 
density profile. The K62 profile is given as: 



VI + (r/r c ) 2 ^l + {n/r c y 



(7) 
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S = k 



1 - 



V 1 + (rt/r c y 



(8) 



where S is the stellar surface-density, r the radius along 
the major axis, r c and rt are the core and tidal radii, 
respectively, and So is the central surface-density. The 
best-fit K62 model, shown in Fig. [51 has a reduced x 2 
value of 1.35, and is described with the following param- 
eters: S = (137 ± 54) arcmhr 2 , r c = (6.21 ± 0.95)' 
and rt = (11.70 ± 0.87)'. Here the uncertainties corre- 
spond to the standard deviation of the fitted parame- 
ters based on the non- linear least square fit. The core 
and tidal radii correspond to (460 ± 75) pc and (860 ± 
86) pc, respectively, a ssuming a distance of 254+*g kpc 
(|Bellazzini et al.l2004T ) to Leo I. These parameters are de- 
pendent on an accurate measurement of the Leo I center. 
It is possible that the central coordinates derived in the 
previous section do not represent the absolute minimum 
of the gravitational potential of Leo I because the stellar 
surface-density may not follow a 2-dimensional Gaussian 
distribution. If we repeat the derivation for determin- 
ing the center of Leo I described in the previous section, 
but define the center to be the position where the stellar 
surface-density has its maximum and fit the K62 model 
to those data, then the core and tidal radii change by 
less than 2%. 

Our best fit K62 parameters are somewhat differ- 
ent from the results reported in the literature (for an 
overview see IH95). While the tidal radius is in agree- 
ment with previous results (IH95 derive r t = (12.6±1.5)' 
and list previous results in their Table 8), our core radius 
is substantially higher. r c is larger by ~ 50% than the 
core radius reported in IH95. This is most likely due to 
the photometry extraction package utilized here, which is 
specialized for crowded field photometry, as well as the 
applied robust completeness corrections, which are the 
largest in the central (i.e. most crowded ) regions. 

We also find the best fit Plummer law (|Plummerlll911f ). 

1 + (r/rp) , and show it in Fig. [9l The re- 
duced x 2 °f the fit is 3.55, and the central surface density 
and characteristic scale are Sq = (149 ±8) arcmin -2 and 
rp = (4.7 ± 0.2)'. The Plummer model is a less good fit 
to the data than the K62 profile. 

In Fig. [5] we also show the best fit Gaussian profile to 
the data; S = Sf exp [-r 2 /(2r 2 )] . The projected stellar 
density is well presented (x 2 = 0.75) by a Gaussian out 
to a radius of ~ 12', with the best fit parameters of 
S^f = (135 ± 7) arcmin^ 2 and r = (3.5 ± 0.1)'. We 
utilize this profile in Sec. [6] to model the mass of Leo I. 

5.3. Tidal tails? 

It is clear from Fig. [9] that the stellar surface density 
falls off less sharply than the K62 model. Such an excess 
of stars may indicate tidal extension; on the other hand 
it may indicate that the (tidally-truncated) K62 model 
is simply a poor fit to the surface-density profile of a 
tidally-undisrupted dSph. Some insight into this issue 
can be gleaned from exploring the 2D density of stars in 
the outer parts of Leo I. 

We display the contour plot of Leo I in Fig.[TOJ This fig- 
ure was constructed using the CMD-selected stars, where 
the background stellar density stated above in £15.21 has 



y p 



been subtracted. Each contour level corresponds to a 
surface-density point in the radial profile. Thus, the 
apparent surface-density excesses at large radii seen in 
Fig. [9] can be directly related to the structures in Fig. 
[T0l For clarity, the contour levels have been divided 
into two groups: those which are well fit by the K62 
model (that is, the points in Fig. [TO] at r < 10'; filled 
contours) and those which deviate significantly from the 
model (r > 10'; solid and dotted contours). The first 
four points beyond r = 10' are represented by solid con- 
tours, and they describe a roughly elliptical shape. That 
is, although the radial profile of Leo I appears slightly 
'inflated' in its outer regions, there is no evidence of a 
distorted structure (such as the 'S'-shaped outer struc- 
ture of the Ursa Mino r dSph; [Marthicz-Dclgad o et al.l 
120011: iPalma et al.ll2003l ). This argues against any strong 
tidal disruption of Leo I. 
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Fig. 10. — The spatial distribution of possible Leo I candidates 
(see color-magnitude filter in the bottom left panel of Fig. [Sjl with 
subtracted field contribution presented in contours. The contour 
levels correspond to surface densities derived in elliptical annuli of 
1.2' width (see text for details). The filled contours are equivalent 
to the surface densities well fit by the K62 profile (first 9 points 
in Fig. [9j . The solid and dotted contours correspond to densities 
poorly fit by the K62 model (the last 5 points in Fig. [9] see text for 
details). The lowest density levels (dash-dotted and dotted lines) 
seem to show asymmetries. The significance of these contours is 
however below 3<r in the background fluctuation of the field star 
population. 

The two lowest contours (corresponding to the two out- 
ermost points in Fig. [§]) are displayed with dotted and 
dashed lines, respectively. A small surface density en- 
hancement is visible to the north-east of Leo I. How- 
ever, this possible excess of Leo I stars is consistent with 
less than 3a in the background fluctuation of the field 
star population. Hence, within the limit of our data set, 
there is no clear evidence for a possible tidal disruption of 
Leo I. This result is not surprising. The relatively large 
luminous mass and Galactocentric distance of Leo I ar- 
gue against its displayi ng st rong signs of ti dal disruption. 
Both fevrd et all (|1994[ ) and lPeeblesI (|1995t ) have modeled 
the orbit of Leo I and argued that the galaxy had at most 
one encounter with a large Local Group galaxy in the 
past. This is consistent with its rather small tidal radius 
and no evide n ce for disruption beyond it. Independently, 
iBowen et al.l (|l997j) showed that there is no evidence for 
tidally disrupted gas in Leo I based on three QSO / AGN 
spectra which they utilized to search for absorption by 
gas within the halo of the dSph. 
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We proceed in our analysis with the conclusion that 
Leo I is not tidally disrupted, at least to a magnitude 
limit of i — 22 and a stellar surface-density of 4xl0 -3 
of the central surface-density. Large-area coverage with 
deeper photometry would be needed to resolve whether 
there are tidal extensions of Leo I at fainter levels. 



1000 



6. 



TOTAL LUMINOSITY, MASS AND MASS-TO-LIGHT 
RATIO OF Leo I 



6.1. Total luminosity 



The CMD-filtered dataset was used to measure the lu- 
minosity function (LF) and total luminosity of Leo I. 
First, we derived LFs for the field and core regions of 
the dSph by counting the number of stars in bins of 0.15 
magnitudes down to the completeness limit of the sur- 
vey, i — 22. In i i5.1l we derived the structural param- 
eters of Leo I, and these were used to define the areas 
from which to draw the core and field populations. The 
core population was taken from within an ellipse with a 
semi-major axis of 12' (approximately the tidal radius), 
while the field population was defined to be that beyond 
a semi- major axis of 30'. Both functions were corrected 
for incompleteness using the estimates listed in Table [TJ 
The Leo I LF was derived to be the difference between 
the core and field functions (after scaling the field func- 
tion to match the area from which the core population 
was drawn). 

Fig. [11] (top panel) shows the completeness-corrected 
and background-subtracted i band luminosity function 
of Leo I. The error-bars are the combination of Pois- 
son uncertainties for the field and core populations, in- 
cluding those taken from the completeness corrections. 
Our magnitude limit of i = 22 excludes part of the red 
clump and the horizontal branch of Leo I. iHeld et all 
(2000) have shown that the morphology of the Leo I 
horizontal branch is remarkably similar to that of the 
intermediate-metallici ty globular cluster M5 (NGC5904; 
ISandquist et alJfl996f ). Hence, to account for the miss- 
ing flux of stars fainter than our i band cutoff we sup- 
pl emented the LF of Leo I with the one for M5 given 
in ISandquist et all (|1996l ). This step required a conver- 
sion of the Leo I luminosity function to the Johnson- 
M organ-Cousin syste m via the empirical relations given 
in lSmith et~aTl (l2002f). Distan ce moduli of 22.02 ± 0.13 
(Leo I: iBellazzmi et all l2004h and 14.41 ± 0.07 (M5; 
ISandquist et al.lll996f ) allowed us to place the LFs on 
the same magnitude scale. Finally, the functions were 
aligned by minimizing the difference in the overlapping 
magnitude range. To estimate the errors of such a proce- 
dure we scaled the two LFs using the overlapping range 
(i) only to J c = 20.55 (corresponding to i = 21) since 
our mean photometric errors beyond this limit are quite 
high, i.e. 0.08 ± 0.02, and; (ii) with no upper or lower 
limits. 

Furthermore, we accounted for the missing flux be- 
yond I c = 27.5 using the completeness c orrected LF 
for NGC6397 given in iPiotto et all (fl99l . We chose 
NGC6397 since it has a metallicity comparable to the 
metallicity of Leo I. We used the dist ance modulus for 
NGC6397 given in IPiotto et ail (|1997l ) and scaled it to 
the LF of M5 which has been previously scaled to the LF 
of Leo I as described above. The resulting supplemented 
LF of Leo I is shown in the bottom panel of Fig. [TTJ 
It is in excellent agreement with the LF of Leo I pre- 
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Fig. 11. — Top panel: Leo I i band luminosity function (LF) cor- 
rected for the contribution of field stars and completeness (see text 
for details). The indicated errors correspond to Poisson and com- 
pleteness uncertainties for field and Leo I stars added in quadra- 
ture. The width of the magnitude bins is 0.15 mag. Bottom 
panel: Johnson-Morgan-Cousins I c band luminosity function of 
Leo I (bold line) supplemented with the LF of the intermcdiate- 
meta llicity globular c l uster M5 (NGC5904; thin dash -dotted line) 
fro m l Sandquist et all Ill996fl and NGC6397 given by IPiotto et all 



(|1997T ) (thin dotted line). The first LF was scaled to the distance 
and number counts (N) of Leo I and the latter to the extrapolated 
LF of Leo I (see text for details). 



sented in lLee et al.l (|1993l ) out to a magnitude of J c ~ 23, 
beyond which their counts are dominated by incomplete- 
ness. By integrating over the luminosity function we cal- 
culated the total I c absolute magnitude of Leo I to be 
Mj c = —12.03 ± 0.14. Assuming a solar absolu te magni- 
tude of M/ c ,0 = 4.14 (jOdenkirchen et al.ll2001h yields the 
total luminosity of Leo I of L Ic = (3.0 ± 0.3)xl0 6 L/ CI0 . 

6.2. Total mass estimate and mass-to-light ratio 

We estimate the central mass-to-light ratio, and on this 
basis the total mass, of Leo I using the Jeans equation 
for an isotropic spherical stellar system (c.f. Binney & 
Tremaine 1987): 



1 d (crip* 



GM(< r) 



(9) 



p* dr r 2 
where r is the radius, <7* is the stellar velocity dispersion, 
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p» is the stellar density. The right hand term in eq. [9] 
represents the gravitational force of the system given in 
terms of the total (i.e. stellar and dark matter) mass, 
M, within a radius r. For a simplified solution of eq. [9] 
we assume that the whole stellar system resides in a dark 
matter halo which dominates the mass. 

To solve eq. [5] analytically for pdmh we need to find a 
simple model that describes the stellar volume-density 
distribution. In H5.2I we have shown that the stellar 
surface-density (that is, the projected stellar volume- 
density) profile is well fitted by a Gaussian function 
(out to a radius of ~12'). Since the projection of a 3- 
dimensional Gaussian distribution is again a Gaussian, 
we adopt p*(r) — p*,o exp(— r 2 /2r 2 ). We also make the 
ansatz that er* w constant, which leads to a solution of 
eq. [S] with a central mass density: 

^ DM '° = (Ir) G% (10) 

iMateo et al.l (|1998bD reported a velocity dispersion for 
the core of Leo I of (8.8 ± 1.3) km s _1 based on 33 red 
giants in the dSph. The average angular separation of 
the observed stars from the center of the galaxy was 1.9' 
with a largest separation of 3.5'. They have also shown 
that there is no significant variation of the velocity dis- 
persion with stellar type or radius. Therefore, the ansatz 
that er* is constant is a reasonable assumption, at least in 
the central part of Leo I. Hence, within a radius of ~3.5' 
around the center of Leo I the dark matter mass density 
is Pdm,o = (0.07 ± 0.02)M©pc~ 3 which corresponds to 
a central M/L ratio of ~ 3 in solar J c band units. Here 
the reported uncertainties are formal statistical uncer- 
tainties. 

The inferred central mass density i s comparable to the 
central den sity of the Forna x dSph (|Mateo et al.l 1199 lh 
and Leo II (|Vogt et al.| [l995) . This is consistent with the 
DM central densities for low-luminosity, gas-rich dwarf 
irregulars (for a review see Mateo et al. 1991). On 
the other hand, it is lower than the central densities 
inferred for other Local Gro up dSp h galaxies (see e.g. 
IMateo et al1l!991l and lMated[l998al for a review). Such 
a result may argue in favor of DM in Local Group dSphs 
being governed by Galactocentric distance; we discuss 
this in more detail in $7] 

With the limited kinematic data-set in hand we can 
only bracket the possible M/L ratio of Leo I (within the 
tidal radius) by investigating two limiting cases: 

(i) We assume that the total mass profile follows the light 
distr ibution. Although t his assumption is n ot very likely 
(e.g. IMateo et al.l Il991h iKlevna et al.ll2001h . it provides 
a good constraint to the minimum mass of Leo I. The 
total mass is then given by: 

/>oo 

M tot = 47r / p maas (r)r 2 dr = (2tt) 3/2 pdm^o (H) 
Jo 

The total mass of Leo I is then (1.7±O.2)xlO 7 M , with 
quoted formal errors. Using the total luminosity of Leo I 
derived in the previous section we infer a mean I c band 
mass-to-light ratio of ~6 in solar units. 

(ii) We assume that the mass-density stays constant to- 
wards large radii, i.e. p mass (r) = Pbm.o- The total mass 
within a given radius is then: 



M (r) = 4tt / pdm.o r' 2 dr' = — p DM .o r 3 (12) 
Jo 

In this case the total mass of Leo I within r = 12' is 
(20 ± 6) x 10 7 M Q and the mean mass-to-light ratio in the 
I c band is ~75 in solar units. 

Hence, for the two limiting descriptions of the distri- 
bution of mass in Leo I we infer an I c band M/L ratio 
of 6 and 75 in solar units, respectively. If we adopt the 
mean V — I c co lor for Leo I of roughly 0.6 (see CMDs in 
iLee et al.lfl993l and lCaputo et al.lll999f) . then the Uband 
M/L ratio is about 5.5 and 65 in solar units, respectively, 
for the two limiting cases. 

7. DISCUSSION 

We have estimated the mass-to-light ratio of Leo I in 
the I c band to be in the range of 6 — 75 in solar units. 
The lower limit was calculated under the assumption 
that mass follows light. The upper limit represents the 
M/L ratio within ~ 12' of a system with a large dark 
matter halo extending far beyond the limiting radius of 
the visible matter, where the dark matter halo density is 
constant within the luminous limiting radius. However, 
given the current evidence concerning dark matter ha- 
los in dwarf galaxies, it is possible to put constraints on 
which of these scenarios is the most likely. 

If the distribution of mass follows the distribution of 
the visible component, then the velocity dispersion of 
Leo I would be expected to fall to zero at the tidal ra- 
dius. Nonetheless, past studies have shown that such a 
behavi our is not typi c al for Local Group dSphs. For ex- 
ample, IKlevna et al.l (|2001h have shown that the radial 
variation of the velocity dispersion of the Draco dSph 
is flat within the uncertainties. Additional evidence for 
a flat velocity dispersion profile is given by IMateo et al.l 
(|1991h for the Fornax dSph. The evidence is similar for 
Leo I. Koch et al. (2007) report that the velocity disper- 
sion profile of Leo I appears constant at increasing radii. 
This suggests that mass does not follow light in Leo I and 
favours option (ii). We infer that the Leo I M/L ratio 
is 3> 5 in the V band. The inferred upper value of the 
M/L ratio in the Uband is ~ 65, yet this represents the 
M/L ratio only within ~ 12'; hence the possibility that 
the true M/L ratio of Leo I is even greater than this is 
not ruled out. 

It has been argued (for example, iKlessen &; Kroupal 
1998) that distortion due to the Galactic tidal field can 
reproduce the appearance of a system dominated by dark 
matter. Under this scenario, the gravitational potential 
of the Galaxy has heated the internal structure of an 
orbiting satellite, thereby leading to an inflated velocity 
dispersion. If we assume a stellar M/L of 3 for Leo I, then 
the velocity dispersion must be inflated by a factor of 4—5 
to produce an apparent M/L of 65. Given that Leo I 
is the farthest known dSph from the Galactic center, it 
would be surprising to find tidal forces strong enough to 
alter the internal kinematics to such a degree. Indeed, 
our analysis found no evidence of tidal disruption. In 
M5.3I we demonstrated that the structure of Leo I follows 
a smooth, elliptical shape to the limiting radius. Thus, 
it seems unlikely that tidal forces are responsible for the 
large measured M/L value, and this result suggests that 
Leo I is strongly dominated by dark matter. 
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8. SUMMARY 

We have presented an automatic photometric pipeline 
especially designed for SDSS crowded-field images. The 
software performs extremely well on crowded SDSS im- 
ages, yielding high-quality photometry with complete- 
ness above 80% down to a magnitude of < 21 and stellar 
density of up to ~ 200 arcmin~ 2 . We have extensively 
tested the detection efficiency and photometric accuracy 
of this software (as a function of both magnitude and 
stellar density). The pipeline was applied to a region of 
~ 3.55D° around the dwarf spheroidal galaxy Leo I in 
the three bands g, ° and i. 

We constructed a filter in colour-magnitude space tar- 
getting the Leo I red giant branch. It removed ~ 80% of 
the foreground contamination, and enhanced the central 
stellar density contrast by a factor of > 4. We find that 
the projected spatial structure of Leo I is ellipsoid-like 
with core and tidal radii (following an empirical King 
model) of (6.21 ±0.95)' and (11.70 ±0.87)', respectively. 
This corresponds to (460 ± 75) pc and (860 ± 86) pc, re- 
spectively. The radial profile deviates slightly from the 
King profile towards the outer regions, however there is 
no evidence for extra-tidal structures (such as tidal tails) 
down to a magnitude of i — 22 and at a stellar surface- 
density of 4xl0~ 3 of the dSph's central density. 

The luminosity of Leo I was measured by integrating 
the observed completeness-corrected and background- 
subtracted flux of stars. We accounted for the missing 
flux to the faintest levels (ic =34) supplementing the lu- 
minosity function with the LFs of the globular clusters 
M5 and NGC6397. The I c band absolute magnitude of 
Leo I was measured to be Mj c = —12.03 ± 0.14 and the 
total luminosity (3.0 ± 0.3) xlO 6 L IetQ . 

We modelled the mass of Leo I using the spherical 
Jeans equation and assumed that the stellar velocity 
dispersion is isotropic and spatially constant. The in- 
ferred central mass density of Leo I is then (0.07 ± 
0.02) Mq pc -3 . Assuming that mass follows light the to- 
tal mass of Leo I is (1.7±0.2) xlO 7 Mq. This value is com- 
parable to the total mass estimates using the standard 
'core' fitting method. Under the assumption that the 
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dark matter halo is extended and its density is spatially 
constant, the total mass of Leo I within the tidal radius 
of the visible component (i.e. 12') is (2 ± 0.6) xlO 8 M Q . 
Hence, the total mass to light ratio for Leo I in the I c 
band and in solar units is in the range of about 6 to 75, 
where the first and latter values correspond to the above 
quoted masses, respectively. In Sec. [7] we argued that 
the mass-follows-light assumption is not reasonable for 
the Leo I system and therefore concluded that the M/L 
ratio must be ^> 6 in I c band solar units and possibly 
> 75 if a constant density DM halo would dominate the 
mass and extend further beyond 12'. 
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